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Abstract 

The NA48 collaboration has measured the amplitude of the CP-conserving component of 
the decay vr+vr^vr'^ relative to 7r"'"7r~7r''. For the characteristic parameter A, the 

values Re A = 0.038 ±0.010 and Im A = —0.013 ±0.007 have been extracted. These values agree 
with earlier measurements and with theoretical predictions from chiral perturbation theory. 

1 Introduction 

The — s> Tr+Tr^Tr*^ decay amplitude is dominated by two angular momentum com- 
ponents, / = and / = 1 (/ is the angular momentum of the neutral pion with respect 
to the system of the two charged pions); higher angular momentum states are suppressed 
because the kaon mass is close to the three pion mass. In this analysis we have measured 
the CP conserving transition to the I = 1 state through its interference with the domi- 
nant —>■ TT^TT^TT^ decay. We neglect any effects from CP violation in mixing and decay, 
which are totally negligeable at our level of sensitivity. 

The amplitudes for the decays of neutral kaons into three pions {Af^ for Kl and 
for Kg) can be parameterized in terms of the Dalitz variables X and Y, which are 
defined as 

- — 2 — ' ^ - — 2 — y^) 

m + m + 

with = {pK — V-kY 1 "^0 = |('57r+ ± s^,- ± s^o), and and being the 4-momenta of 
the kaon and the pion respectively. The Dalitz variable X is a measure of the difference 
of the energies of the two charged pions in the kaon's rest system while F is a measure of 
the energy of the neutral pion in the kaon's rest system. 

The / = 1 and the / = components of the decay K\ 7r"'"7r~7r° can be separated 
by the fact that the amplitude of the / = 1 process is antisymmetric in X while the 
/ = process is symmetric in X. Therefore the / = 1 contribution can be extracted by 
separately integrating over the regions X > and X < and subtracting the results from 
each other. We consider the distribution 



V(t) = 3^ (2) 



where N^'^^it) [X^^°(t)] is the number of decays of neutral kaons into tt+tt^tt'^ at time 
t with a value of the Dalitz variable X larger [smaller] than zero mi21- 

In the NA48 set-up, where Kg and K^ are produced in equal amounts at a fixed 
target, the interference of the / = 1 component of Kg tt+tt^tt*^ with the / = component 
of the dominant K^ — ^ vr+vr^Tr^ decay can be observed. It can be described by a complex 
parameter j2] 

^ ^ IZo dY r dX Al ^-^'=°^ (X, y )4-^'=^^ (X, Y) 
fZ^dYf^-dX\Al-^'='\x,YW 
which has been extracted by fitting the distribution defined above: 

_ 2D{E)[Re{X)cos{Amt) - /m(A)sm(Amt)]e"^^^"^^^ 

V [t) ~ — ^- . [4:) 

e 

Am is the mass difference between K^ and Kg, and Tg are the respective lifetimes, 
and the energy-dependent "dilution" D{E) is the difference in the relative abundances of 
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and K at production for a kaon energy of E: 



K° + K 

2 Experimental setup 

The data presented here were taken in a high-intensity neutral beam at the CERN 
Super Proton Synchrotron during 89 days in 2002. Protons of 400 GeV with an average 
intensity of 5 x lO^'' per 4.8 s spill (within a cycle of 16.2 s) impinged on a beryllium 
target, where secondary particles were produced. Charged particles were removed by a 
sweeping magnet and a 5.1 m long collimator, which selected a beam of neutral particles 
at an angle of 4.2 mrad with respect to the proton beam. 

After the collimator, the neutral beam and the decay products propagated in an 
89 m long vacuum tank, where most of the short-lived kaons (K^) and neutral hyperons 
(A*^, 'EP) as well as a small fraction of the long-lived neutral kaons {K^) decayed. Undecayed 
particles {K^, neutrons) continued to a beam dump via a vacuum pipe that passed through 
all of the downstream detectors. 

Downstream of the vacuum tank there was a magnetic spectrometer consisting of 
four drift chambers inside a helium-filled tank and a magnet with a horizontal transverse 
momentum kick of 265 MeV/c. The spatial resolution was 120 /im per view, and the 
momentum resolution was 6p/p = (0.48 © 0.015 p)%, with momentum p in GeV/c. 

After the spectrometer there was a scintillator hodoscope, for the accurate timing 
of charged particles (yielding a time resolution on single tracks of 250 ps), and a liquid- 
krypton electromagnetic calorimeter, for the measurement of photons and electrons, with 
an energy resolution of a{E)/E = {3.2/VE © 9/E © 0.42)% {E in GeV). These were 
followed by a hadron calorimeter and a muon detector. A more detailed description of the 
detector has been given elsewhere |3]. 

3 The trigger 

Decays into vr+vr^Tr^ made up only a small part of the event rate in the detector. This 
was due to the small expected branching ratio for the Kg decay and to the long lifetime 
of the K^. To obtain a statistically significant sample of tt+tt^tt'^ decays, a very selective 
trigger was used to suppress the high two-particle background, which came mostly from 
Kg TTTT and from A*' — >■ p7r~ . 

When the level-1 trigger indicated an event with at least two charged tracks and a 
minimum energy deposition of 30 GeV in the calorimeters, a fast online processor calcu- 
lated the momenta of the charged tracks, and hence the invariant mass of the decaying 
neutral particle under the assumptions of Kg tt+tt^, A° pvr^, and A — pvT"*". Events 
were accepted by the trigger if none of the calculated masses matched the physical mass 
of the particle in question (difference of squares of the calculated mass and the nominal 
mass [5. larger than two percent). Additional cuts were applied to the relative momenta of 
the two charged particles {piarger I V smaller < 3.5, to further suppress A° decays) and to the 
separation of the two tracks at the entry of the spectrometer (> 5 cm in drift chamber 1, 
to remove photon conversions in the upstream Kevlar window) . To allow a measurement of 
the trigger efficiency, a downscaled sample of minimum-bias events (1/35) was recorded in 
parallel, where only the level-1 trigger condition was fulfilled. In the course of data-taking, 
the high voltage in the drift chambers was occasionally adjusted to maintain stable oper- 
ating conditions. As a result, the efficiency of the level-2 trigger varied between (75±4)% 
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and (86±4)%. By comparing the analysis of different time periods, it was cliecked that 
this variation did not introduce any spurious asymmetries. 

4 Data analysis 

For the offline analysis, events were selected with two tracks of opposite charge from 
the same vertex (distance of closest approach less than 3 cm), and at least two electro- 
magnetic clusters (within ±20 ns). For each pair of clusters, the vr^-mass was calculated 
assuming that the pion decayed at the position of the charged vertex. If there were more 
than two clusters, the cluster pair yielding the best 7r°-mass was retained. 

In order to reduce acceptance effects near the edges of the detectors, each charged 
track was required to cross the drift chambers at a minimum distance of 12 cm and a 
maximum distance of 110 cm from the beam axis. Likewise, each photon cluster was 
required to be more than 15 cm from the beam axis and more than 10 cm from the outer 
edge of the calorimeter. 

In addition, the photon clusters had to be well separated from tracks (>15 cm) to 
avoid mismeasurements due to energy sharing. Tracks were identified as charged pions 
rather than electrons by demanding that their electromagnetic energy deposition was 
less than 90 percent of the track momentum. Cuts were applied on the charged vertex 
(less than 3 cm lateral distance from the beam line, and between 800 and 5400 cm 
from the target), and on the masses reconstructed for — >■ 7r"'"7r~7r° and 7r° 77 
(m^?^ ± 10.5 MeV and m^i^^ ± 7.8 MeV 0). These cuts also excluded K^g ^ 
decays. The reconstructed kaon energy was required to be between 30 and 166 GeV. A° 
hyperons were further suppressed by cutting harder than in the trigger on the masses for 
the hypotheses A° — pn~ and A pvT"*" {m^^'^±21 MeV |SI) and on the momentum ratio 
of the two charged particles {piarger/v smaller < 2.9). No A° hyperons survive after these 
two cuts. The mass cuts were optimized to maintain a high efficiency without introducing 
substantial background. A total of 19 million events passed all these cuts. 




200 



cm GeV 
z vertex position Kaon energy 

Figure 1: Position of the decay vertex along the beam line (left, measured from target 
position) and reconstructed energy (right) for the data sample of 19 million events. 

Fig. shows the distributions of the position of decay vertices along the beam line 
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(z vertex) and the reconstructed energy after applying the analysis cuts. 




Figure 2: The invariant mass. The vertical lines indicate the K'^ mass cut used in the 
analysis. 

Fig. 121 shows the distribution of the invariant mass after cuts. 

Applying these cuts and correcting for acceptance and trigger efficiency as described 
below, the experimental distribution V(t) defined in (j21) was obtained. The values of Re A 
and Im A were extracted by fitting equation (0]) to the experimental data. Due to the 
energy dependence of the dilution D{E), the fit was made in bins of kaon energy (see 
Fig. El), with 8 equal-sized energy bins ranging from 30 to 166 GeV. The values for the 
external parameters Am, and ts were taken from |5'. For the dilution D{E), we used 
a quadratic fit to the values measured in ^ (energy E given in GeV): 

6 XO 1 T2 

DmeasurediE) = -0.167 + j^E - -^E^ (6) 

The errors for the individual points measured in |ni vary between 5 percent at the lowest 
kaon energies (75 GeV) and 25 percent at the highest energies (165 GeV). We assume the 
uncertainties at the different energies to be uncorrelated. To take into account the differ- 
ence in proton energy and production angle between ^ and our own measurement, we 
used measurements with charged kaons J7], from which we derived an empirical correction 
factor (energy E given in GeV): 

^ ^ 103 105 10^ 1010 ^ ^ 

To fit equation (HI), we then use 

D{E) = Dmeasured{E) X Dcorr{E) (8) 

Based on the uncertainties in and the overall relative uncertainty in D{E) is 
estimated at about 15 percent. 
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The aim of the experiment was to measure a small effect that is antisymmetric in 
the Dalitz variable X (the interference of and K^^ in the decay into tt+tt^tt^) over a 
large background that is symmetric in X (from the CP-conserving decay — > 7r"*"7r~7r*'). 
Therefore, care had to be taken to exclude any detector effects that could artificially 
introduce such an asymmetry in the decay. 

The sign of X equals the sign of the momentum difference of the 7r+ and the vr" 
in the kaon rest system. The pion momenta in the kaon rest system are correlated to 
their momenta in the laboratory system. Any momentum-dependent difference in the 
acceptance or the trigger efficiency for positive and negative pions may therefore introduce 
a possible bias. Due to the field of the spectrometer magnet, positive and negative particles 
populated the detector in a different way. In order to compensate such effects, the field 
of the spectrometer magnet was reversed weekly during the data taking. 

An artificial asymmetry can be caused by the detector acceptance only if this accep- 
tance is asymmetric in X. However, if there is a real, physical asymmetry, its value may 
be distorted even by a symmetric X-dependence of the acceptance. The NA48 detector 
simulation (based on the GEANT3 package [Hj) showed a symmetric decrease in the ac- 
ceptance for high Using the simulation, the acceptance was calculated as a function 
of kaon decay time and kaon energy. The effect was then corrected by weighting each 
event in the data in inverse proportion to the calculated acceptance for the corresponding 
value of kaon decay time and kaon energy. This procedure resulted an increase of 
10 ± 5 percent in the extracted values of Re A and Im A. 

In the drift chambers, signal attenuation near the ends of some wires created a 
left-right acceptance asymmetry on the detector periphery, especially for strongly devi- 
ated low-momentum particles. This was compensated in the offline reconstruction by the 
redundancy of chamber planes, and there was no visible effect on the minimum-bias data 
sample. However the data collected with the fast level-2 trigger algorithm, which did not 
use the full redundancy, showed an asymmetry in X at all kaon lifetimes, which was 
strongest (2 percent) at the lowest kaon energies (30 GeV). This asymmetry was reversed 
when the magnetic field of the spectrometer magnet was inverted. Using the X-symmetric 
amplitude of Kl — >• 7r"'"7r~7r°, the detector and trigger simulation accurately reproduced 
this behavior. 

The Monte-Carlo distributions from the simulation were used to correct the data. 
The correction resulted in a shift of about 3 x 10~^ in both Re A and Im A, with opposite 
sign for the two field orientations of the spectrometer magnet. The simulation program 
had been developed, extensively used and validated during previous studies by the collab- 
oration lUIH]. After acceptance and trigger corrections, the individual results for the two 
field orientations differ by about 1.4 a. To exclude possible additional effects the Monte- 
Carlo corrected data samples with positive and negative fields in the spectrometer were 
subsequently averaged (taking the arithmetic mean of the observed distributions V{t) in 
formula which is equivalent to normalizing the data taken in the two field orientations 
to the same number of events and summing them before applying (0)). Slightly different 
results were obtained when using either only the Monte-Carlo correction (without av- 
eraging over the two magnetic field orientations), or only averaging over the magnetic 
fields (without the Monte-Carlo correction). We use half the difference between these two 
last approaches as a contribution to our systematic error (±0.004 in Re A and ±0.003 in 
Im A). 
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Figure 3: The experimental distribution V{t) (points with error bars) and the function 
obtained by a fit (smooth curve) in bins of time and energy. The P^i' degree of 
freedom is 169.35/162=1.045. The fit was made simuhaneously over all energy bins, but 
is mainly constrained by the intermediate energy bins. The data were compensated for 
acceptance and trigger efficiency by applying Monte-Carlo corrections and by averaging 
over samples with opposite magnetic field in the spectrometer. 
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Another cross-check was to analyse the data with a tight outer radius cut of 50 cm 
in the drift chambers and otherwise follow the above procedure. This was prompted by 
the fact that the observed drift chamber inefficiencies were highest at large radii. The 
results in Re A and Im A differed by about 0.001 while the statistical errors increased by 
about the same amount. 

The other major contribution to the systematic uncertainty comes from the /K 
dilution D{E) defined in formula (0) and leads to an additional uncertainty of ±0.005 
in Re A and of ±0.003 in Im A. By varying the experimental cuts quoted above within 
reasonable bounds and the external parameters from ^ within their errors it was seen 
that all other possible sources of systematic uncertainties were negligible. 

^.02 
0.015 
0.01 
0.005 




2 4 6 8 10 
Kaon decoy time t/tj 




Figure 4: The experimental distribution V{t) (points with error bars, data summed over 
all kaon energies) and the function obtained from the fitted values for Re A and Im A for 
a mean kaon energy of 87 GeV. 



5 Results 

The following values have been obtained for the real and imaginary parts of A: 

Re X = ±0.038 ± 0.008,^* ± O.OOGsyst (9) 
Im \ = -0.013 ± 0.005,tat ± 0.004,y,t (10) 

Adding both errors in quadrature, 

i?e A = ±0.038 ± 0.010 (11) 
Jm A = -0.013 ±0.007 (12) 

These values can be compared with the theoretical numbers of Re A=±0.031 and 
Im A=-0.006 obtained in ^] by using results from Chiral Perturbation Theory 

Fig. m shows the distribution V{t) defined in Q derived from the data, and the 
curve obtained from with the fitted values of Re A and Im A (central curve, plotted 
for a mean kaon energy of 87 GeV; the top and bottom curves result from varying the 



7 



dilution within errors). These values were not obtained by fitting this histogram directly, 
but by fitting all data in bins of energy. This is important because of the dependence of 
the dilution D{E) on the kaon energy. 

The real and imaginary parts of A are strongly correlated, as can be seen from Fig. 
which shows the measured values and contours of equal the complex plane. Taking 
into account only the statistical errors, the correlation coefficient between Re A and Im A 
is 0.66. 




0.02 0.04 0.06 



Re A 

Figure 5: The fitted values of Re A and Im A in the complex plane with the contours for a 
of 1 to 10 above the minimum x^- The shorter error bars and the contours refer only to 
the statistical uncertainty while the longer error bars correspond to the total uncertainty. 

From the real part of A, the branching ratio for the CP-conserving component of the 
decay — » 7r'*'7r~7r° can be obtained. Using the same approximation and the external 
parameters as in |3j, we get: 

BR{K'I tt+tt-ttO) = (4.7 tf^jistat) tli{syst)) x 10"^ 

The results from this experiment agree with chiral perturbation theory and 
with two other measurements with comparable errors E|. It is worth noting that 
experiment used tagged kaons, so that its result is independent of any measurements 
or calculations of the dilution D{E) while ^U] was carried out at higher kaon energies 
and therefore higher dilution than our measurement. It should also be noted that all 
three experiments obtain a non-zero negative value for the imaginary part of A (which is 
predicted to be small). 
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